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1. Introduction 



Since tlie pioneering work of de Vega and Woynarovich for the construction of models 
with alternating spins quite a lot interesting generahzations have been presented 0, 
Otherwise, not so much results were obtained concerning the physical structure of the 
models, e.g. the low temperature behaviour of heat capacity and magnetic susceptibility. 
Even the structure of the ground state in the framework of the Bethe ansatz is not fully 
known for the original model. 

In this paper we therefore continue our investigation of the XXZ(|, 1) model with 
strictly alternating spins started in paper which will be lower referred to as paper I. 

In section 3 the thermodynamic Bethe ansatz (TBA) is analysed for zero 
temperature in different regions of coupling constants. 

Section 4 deals with the conformally invariant model, where the low temperature 
behaviour can be determined analytically. 

In section 5 we derive some new results for the isotropic case XXX(|, 1). 

Our conclusions are contained in section 6. 



2. Definition of the model 

We send the reader to papers I and for the basics of the model. We will follow below 
the definitions and notations of paper I. 

Our Hamiltonian of a spin chain of length 2N is given by 

n{-i) = cH{-i) + cH{-i) - HS\ (2.1) 

with the two real coupling constants c and c. Anisotropy parameter 7 is limited 
< 7 < 7r/2. 

For convenience we repeat the Bethe ansatz equations (BAE) and the magnon 
energies: 

/ sinh(A^- + q) sinh(Aj + ^ |t sinh(Aj - + n) ^ 

ysinh(Aj — q) sinh(Aj — ^7) j ^J--^ sinh(Aj — A^ — ^7) ' 



E = cE + cE- -M]H, 



M 



2 sin 7 



^ cosh 2\j — cos 7 

~ 2sin27 
E = 



^ cosh 2\j — cos 27 



(2.3) 
(2.4) 
(2.5) 
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3. Thermodynamic Bethe ansatz and the ground state for different signs 
of the couphng constants 

In section 3 of paper I the TBA was considered for special values of 7 = 7r//i, integer 
and /X > 3. We also argued, that the ground state structure is uniform in our whole 
7-region, while we expect possible changes for the excitations at the 7-points above. We 
therefore use the results of I for the possible appearance of strings in the ground state 
according to the different regions of couplings. 

For completeness we quote in all cases the TBA, equations (3.19) of paper I. We 
found it more convenient to use A-space instead of Fourier transformation, which one 
easily can derive from our equations below. We then recall 

Pf\ ^ 2 sin 727 

/ A,n,±l =± — rTTVT • 3.1 

cosh 2 A =F cos wy 

For shortness we drop the magnetic field in the TBA. Afterwards it can be added without 
any problem. 

Now we analyse the zero temperature TBA in the various regions of signs for c and 

c. 

(i) c > 0, c > 



e+(A) = -c/'(A,l,l)-c/'(A,2,l) 

/'(A, 2,1) 



5{X) 



27r 



r(A,i,i)+m3,i) 

27r 



(3.2) 



4(A) = -c/'(A, 2, 1) - g[/'(A, 1, 1) + /'(A, 3, 1)] 
7'(A,1,1) + /'(A,3,1)' 



27r 



5(A) + 



2f(A,2,l) + f(A,4,l) 
2% 



* £2 • 



(3.3) 



The solution has been already given in 0, where also the excitations have been found, 
(ii) c > 0, c < 

We expect (1,+) and (1,-) strings. 

e+(A) = -c/'(A,l,l)-cf(A,2,l) 

/'(A, 2,1)" 



5(A) + 



27r 



/'(A, 2,-1) 



27r 



1) 



(3.4) 



eli(A) = -cr(A,l,-l)-c/'(A,2,-l) 



[/'(A, 2, -1)1 




2n 


* — 



5(A) 



f(A,2,l) 
2n 



* e 



(3.5) 
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At the first hand, one might expect that the solution is given when both strings are 
distributed with infinite Fermi radius. We have determined this state and calculated its 
energy. But it is not the ground state. The same applies to the state with only (1,+) 
strings. That can be seen already superficially after obtaining Sz ^ ^ for it. 

The situation changes when only (1,-) strings are considered. This is due to the 
fact, that the two last terms in equation (3.4) are definitely non-negative while this is 
not the case in equation (3.5), where the term after the (5-function spoils the argument. 

Equation (3.5) for ej^(A) = has been already solved in I. 

e-i(A) 



TT — 7 cosh(7rA/(7r — 7)) 

Attc cos(7r7/2(7r — 7)) cosh(7rA/(7r — 7)) 



+ 



77 — 7 cosh(27rA/ (vr — 7)) + cos(7r7/ (vr — 7)) 
Introducing the function (7 (A, a) 

4:71 cos(7ra/2(7r — 7)) cosh(7rA/(7r — 7)) 



g{X,a) 



vr — 7 cosh(27rA/(7r — 7)) + cos(7ra/ (vr — 7)) 
the solution of equation (3.4) can be written as 

e+(A) = ~cg{X, 7r/2 - 7) - 5g{X, 7r/2 - 37/2). 

For consistency it is necessary to have 

eli(A) < and e+(A) > 0. 

Both conditions specify the region of c and c where our solution is valid. 
We start with ei(A). 



(3.6) 



(3.7) 



(3.8) 



(3.9) 



ei(0) 



27r 



TT 



7 



+ 



cos(7r(7r — 27)72(71 — 7)) cos(7r(7r — 37)/2(7r — 7)) 
00 one has 



> 0. 



Considering the asymptotics for A 

2tt \_ 7r(7r — 27) _ 71(11 — 37) 

C cos ^ h C COS 



vr 



7 



2(vr-7) 



2(vr-7) 



> 0. 



(3.10) 



(3.11) 



We now assume that the two necessary conditions ( |3.10| ) and (|3.11| ) are also sufficient 
to fulfill the second part of ( p.9|) . 

The smaller one of the ratios of the two cosine-functions is then the upper limit of 



c/\c\. Hence after elementary recasting 

1 



|c| 2 cos(7r7/2(7r — 7)) 



< 7 < 



27r 
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We treat e_i(A) in the same way obtaining 

7r7 



< 2 cos 



2(vr-7)' 



(3.13) 



Now it is not difficult to show that condition ( |3.13| ) is fulfilled, when (3.12) holds. 

Therefore, our solution, a sea of (1,-) strings with infinite Fermi zone, is the ground 
state configuration as long as the inequalities (3.12) hold. In the (c, c)-plane this is an 
open triangle formed by the negative c-axis and the straight line given by relation (3.12) 
when the equality holds (s. figure 1). For 7 — (isotropic case, s. section 5) this is 
c/\c\ = |. For increasing 7 the region first enlarges until 7 = 27r/5 and then shrinks and 
approaches the c-axis when 7 7r/2. 

Above that line we expect still (1,-) strings but together with (1,+) strings. So 
moving counter-clockwise from the positive c-axis towards that line the Fermi radius of 
the strings with positive parity shrinks from infinity to zero, while the radius for the 
strings with negative parity is infinite, as can easily be seen from equations I (3.17), 
which implies in the case H = 0,that its energy function does not change sign and is 
therefore strictly non-positive in the limit T ^ 0. 

It is remarkable, that a finite Fermi zone occurs without the presence of a magnetic 
field. Apparently the second coupling plays the role of an external field. 

(iii) c < 0, c > 

We expect (2,+) and (1,-) strings. 



-(A) = -c/'(A,2,l) 



+ 



5(A) + 

mi 



c[/'(A,l,l) + /'(A,3,l)] 
2/'(A,2,l) + /'(A,4,l) 



27r 
+ /'(A,3, 



2n 



* e 



15 



(3.14) 



eli(A) = -cr(A,l,-l)-c/'(A,2,-l) 

7'(A,1,-1) + /'(A,3,-1) 



271 



* en 



5(A) 



f (A, 2,-1) 

27T 



* e 



(3.15) 



We have found that qualitatively the same arguments apply as in the case (ii) before. 
Thus, we consider first only (1,-) strings with infinite Fermi radius. Now it is necessary 
to assure e^(A) > in addition to the first condition of ( |3.9| ). Instead of condition (|3.13| ) 
it gives now 



|C| ^ 2 

c ~ cos(7r7/2(7r — 7)) ■ 



(3.16) 



which guarantees e_i(A) < 0. When calculating one has to be careful when the 
Fourier transformation of /'(A, 3,-1) is to be taken. It vanishes for 7 = 7r/3 and 
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changes the sign after that point had been passed. Finally one obtains 

e+(A) = -cg{\ Tx/2 - 87/2) - cg{\ 7r/2 - 7) - cg{\ 7r/2 - 27), < 7 < vr/3, 

e^(A) =0, 7r/3 < 7 < tx/2. (3.17) 

There is no contradiction with formulae I (3.24), which gives two different values for 
7 < 7r/3 and 7 = 7r/3, while larger 7- values were not considered there. 

Let us first consider < 7 < 7r/3. Then from equation ( |3.16|) we have the two 
conditions 

27r / c 



vr — 7 Ycos(7r(7r — 37)/2(7r — 7)) 

c 



+ ^, r^^^TT ^ > 0, 



cos('7r('7r — 27)/2('7r — 7)) cos('7r(7r — 47)/2(7r — 7)) 



2-K ( _ 7r(7r — 37) 
c cos 



vr — 7 Y 2{ti — 7) 



7r(7r-27) , 7r(7r-47)\ 

+CCOS ; — + CCOS ; — > 0. 3.18 

2(vr-7) 2(7r-7)y- ^ > 

Straightforward calculation gives 

IcI 7r7 

— > 2 cos — and 

C 2(77 — 7) 

|c| ^ 8cos^(7r7/2(7r-7)) 



c 



4 cos^ (777/2(77 — 7)) — 1 



The upper term of the RHS is always smaller than the RHS of ( |3.16| ). Hence we have 
to find the maximum of the two RHS of formula (3.19) and (|3.16| ). In our 7-region the 
second inequality of (3.19) is the most restrictive one. Putting things together we find 
for the region with (1,-) strings only 

IcI 8 cos'^(7r7/2(7r — 7)) 7r 

-r > y_J±_A ^Z^^ < 7 < -, 

c 4 cos^(7r7/2(7r — 7)) — 1 3 

|c| 2 7r 7r , , 

^ > r, o ^ ^ < 9 • (3-20) 

c cos 777/2(7? — 7j 3 2 

In the (c, c)-plane that is an open triangle formed by the negative c-axis and the straight 
line given by relation (3.20) when the equality holds (s. figure 1). For 7^0 (isotropic 
case, s. section 5) this is \c\/c = |. For rising 7 the region shrinks and approaches the 
c-axis when 7 — > 7r/2. 

Above that region we expect (1,-) strings together with (2,+) strings the latter with 
finite Fermi radius. The picture resembles region (ii) treated before. 

(iv) c < 0, c < 



7 



et^i\) = -cf'{\,l,-l) - df'i\,2,-l) 



6{X) + 



*eZi. (3.21) 



Here the vacuum is formed by (1,-) strings only. 

f( A,2,-: 

This region was studied in I where also the excitations have been found. 

(v) c = 0, c > 

(vi) c > 0, c = 

We add nothing new to both cases considered earlier in and . 

Now we can summarize our results about the the ground state structure for different 
values of coupling constants. There are four regions and two singular lines (v) and (vi). 
In the two regions with equal signs (which contain the line c = c) the ground state is 
independent of the values of c and c. Here also the Fermi radii are infinite. There is no 
mass gap in the excitation spectrum. 

In the two other regions infinite and finite Fermi radii occur and the concrete 
structure of the ground state depends on the ratio c/c. Nevertheless, we expect them 
to be gapless, too. 

The picture is not fully symmetric, because region (i) is separated from all others 
by highly degenerate ground state on both lines. This is connected with the fact, that 
one sort of strings has to disappaear at once. 

Finally, the model shows an antiferromagnetic behaviour everywhere (for vanishing 
magnetic field) as long as 7 > 0. The isotropic case is considered in section 5. 



4. Calculation of the low temperature behaviour in the case c = c 

In this section we calculate the low temperature heat capacity and magnetic 
susceptibility for vanishing magnetic field in the case c = c. 

We therefore go back to equations (3.10)-(3.13) of paper I where T is considered 
to be small but finite. Instead of paper where the free energy was calculated, we 
use a method due to Wiegmann [0, which for our purpose was used by Babujian and 
Tsvelick ^ to obtain the results for the XXZ{S) model from entropy and polarization. 
To explore this method it is necessary to ensure 7 < 7r/3, 7 = 7r//i; /i integer. 

We will present in some detail the case c = c = c < while for c > we mention 
only the necessary changes and the final results. 

For c < we have 

ej > 0, j = l.../i-l, 

and 

for T ^ if j = 1 . . . /i - 1, 
^ for T ^ 0. 
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(1,-) 


C 

(1,+) 




(2,+) 




\ / 






(1.-) 


\ (1,-) 


Figure 1. The phase structure of the XXZ{^, 1) model for 7 = 



7r/3. The ground 

state strings are indicated for the four different sectores. The arrows symbohze the 
decreasing Fermi radii of the corresponding strings. Broken lines are coordinate axes. 
Where axes are drawn solid they coincide with sector borders. 



Our aim is now to recast BAE and TBAE in a form where energies and densities are 
given through their zero temperature hmits e^'^-' and ej^'* and values vanishing for T ^ 0, 
i. 6. the energy functions ln(l + exp(— ej/T)), j — 1 . . . /i — 1 and ln(l + exp(e^/T)). 

The main step is the multiphcation by which after some algebra leads to the 
systems 



e,, = —^e 



f -E^3^*(-ir^'^2^1n(l + exp(-|)) 

-ii'*(-ir(^)Tln(l + exp(^)), 



= Sj.^i^ef + J: Bj, * (-l)'-('^)rin (l + exp 
^ k=i ^ 

-ir*(-ir('^)rin(l + exp(^)) (4.1) 



and 



^Ji-l 



k=l 



-{-iyO){p^ + p^) = —ef + Bjk * Pk + Qj * P„ 



(4.2) 



k=l 



where we have introduced 



K{\) = -T,,*A-'^{\) 

Q.(A) = -r,,*A;;(A) 

BjkiX) = (Tjk + Tj^ * A-'^ * T,k){X). (4.3) 

Now we want to perform a shift of the A- variable in the functions ej{X) in the following 
way: 



if T 
VM) = ( A + aln 



27r c 



(4.4) 



where the constant a will be determined yet. We choose it in a way that (pj{\) for A — > oo 
has a finite hmit if T — > 0. Prom I (3.24) we can see, that for all ^f\X) ~ exp(7rA/(7r— 7)) 
if A — > 00. Therefore a — — (tt — 7)77. 

After that shift the system (4.1) is rewritten as 

V,-^-12Qk* i-^r^'^ ln(l + exp i-^,)) - K * (-l)^('^) ln(l + exp {^,)) 

k=l 



1 / ttA 

H exp 



TT — 7 



TT — 7 



1 + 2 COS 



TT 



2{p-l) 



= Sm-i^ + E Bjk * (-l)'-('^ ln(l + exp i-^,)) - K * (-l)^(^) ln(l + exp {^,)) 



k=l 



TT — 7 



■ exp 



ttA 



TT — 7^ 



X 



X 



COS — — h COS — — + COS — —. -r 

z(/i — 1) 



2{p-l) 2(//-l) 

The system (4.2) is treated analogously. 

After differentiating equation (4.5) one obtains the important relations 

p-l, T 



(4.5) 



Pi A In ■ 



p 27r|c| 



ior j — 1 . . . p. 
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These relations are necessary to make the appropriate substitutions of variables in 
the integrals for S and S^. No such relations are expected as soon as c 7^ c. 

The starting point for the heat capacity calculation is the expression for the entropy 



S 
N 



E 



Pj In 1 



Pi 

Pi. 



+ pj In 1 



Pi 

Pj. 



(4.7) 



Using sjTumetry and Pj/ Pj = e^^''^ we have 

2 E rfA [pj In (1 + e'^/'^) + pj In (l + e"^^/^ 



S 
N 



(4. 



For T ^ the main contribution to the integral comes from A ^ 1. 

After performing the shift (|4.4|) and using relations (4.6) the entropy becomes 



S_ 

N 



7r — 7 



E(-i) 



00 



T 
2^ 



dX 



d_ 
dX 



ln(l + 



In (1 + e 



d 



+— In (1 + In (1 + e" 



(4.9) 



We are interested only in the leading order for vanishing temperature. Therefore we can 
substitute the lower limit by —00 (both integrals converge). 

We will see in a moment that the remaining integral is even independent of T. 

Now it is straightforward to change the variable in the way 

1 



X 



(4.10) 



1 + e^^ 

for every integral in the sum (Note the change in the definition of the function / in 
equation (3.18) of paper I). 
In final form 



N 



vr — 7 



/(<) 



dx 



Inx ^ ln(l — x) 



1 — X 



X 



with Lff = Lpj(±oo). 

The integral is given by the function 7(0, b) already introduced in 
also the necessary special values have been taken. 



(4.11) 



from where 



7(a,6) 



dx 



\nx 1 , . 
+ - In 1 - X 

1 — X X 



(4.12) 



Following the standard procedure it is more convenient to use another form of the 
TBAE to determine ipf. They differ from those of the XXZ(S) model in the terms 
with the coupling constants only. It is the system (3.17) from paper I. 
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After the shift it takes the form 

ipi{\) = -s* ln/(v92)(A) + |c| exp 



ttA 



TT — 7 



^■(X) = -s* ln[/((^j+i)/((^j_i)](A) + |c| exp 



TT - 7^ 



2T 



s * ln/(v9^_2)(A) 



2T 



+ s*ln/(v?^_2)(A). 



For A — s> — oo the inhomogeneous terms generate a solution of the form 
(fj = +00, j = l...n-l, 

which imphes 

/(^7) = o, j = i.../i-i, 

For A — CX3 the free terms can be neglected, and thus the solution is given in 

sinh(i7/2T) 



For H ^ then 



</^+_i = ^ + ln 



smh{H{j + 1)/2T) 

sinh(iJ(/i - 1)/2T) 



2T 



In 



sinh(/7/2T) 
sinh(ij'(/x - 1)/2T) 
sinh(/f/2r) 



1 



/«-i) = -, /«) = i--. 



(4.13) 



(4.14) 



(4.15) 



(4.16) 



(4.17) 



We mention that the above solution does not depend on the sign of the coupling constant. 
The consequences of that fact will be considered below. 
Now we can calculate relation (4.11): 



s 


(vr- 


7)T 


N ~ 


7r2 


c 




(vr- 


7)T 






c 




(vr- 


7)T 




^2 


c 




(tt- 


7)T 



V-2 

E7 



2' 







7 



-,0) -7(l--,l 



+ 1 



3|c| 



(4.18) 
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Finally, for the heat capacity per site (of a chain with 2N sites) 

7r — 7^ 



c 



6|c| 



-T. 



(4.19) 



The polarization is obtained in the same way as in starting with the basic formula 

^ ' ' ' (4.20) 



r /"oo POO 

/ dXp,_,{X)- dXp,{X) 



N 2 

which for our model takes the same form. Using symmetry and the shift gives 



^ = 2^ 
N 2 



/i 



rfAp^_i(A) - 
° (7r-7)T d 



dXp^{X) 



27r2|c| 



\n{l + e'^^-')dX + - / 

J — ( 



(tt - 7)T (9 



27r2|c| aA 



ln(l + 6-'^") dX^ 



/i(7r-7)T 
27r2|c| 



In 



1 + exp(-v?+) 
1 + exp(v9 



In 



1 + exp(-y^ 
1 + exp{ip 



(4.21) 



For the first term we use relations (4.16) with H ^ T. 



_ Hp 
~ 2T 
Hp 



+ In 



2T 



-In 



and hence 



1 + exp(v9+_i) = exp 
1. 



exp{H{p - 1)/2T) 
exp{H/2T) 
exp{H{p-l)/2Ty 
exp{H/2T) 

Hip-1)\ 



H{p-l) 



T 



H 
Y' 



(4.22) 



T 



(4.23) 



1 + exp(-v9;; 

In the second term term (4.14) must be completed by corrections containing the leading 
term in H. 



In 



" 1 + exp(-y?^; 
1 + exp((^-„i;j 
Putting things together 

(tt - ^)H 



H^i 



27r2|c| 

and finally for the susceptibility 

P — 1 TT — 7 



X 



47r|c| 



47r7|c| 



(4.24) 



(4.25) 



(4.26) 



Now we perform the calculation for c = c = c > 0. We will list only the necessary 
changes in the calculation before, induced by the sign of the coupling constant. For 
c > we have 



ei,e2 < 0, 
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The shift in equation ( [4 .41 ) must now be taken as a = — 7r/7 according to the asymptotics 
of the e^^''. Instead of equation (4.6) we have then 



pj I A — /i In 



I A — In 



T 

27rc 
T 

27rc 



Note the change of the overall sign. 
Consequently, (4.11) is modified 



S 
N 



1 



dx 



Inx ^ ln(l — a;) 



1 — X 



X 



(4.27) 



(4.28) 



The change in the system (4.13) is obviously the replacement of |c| by — c. As already 
mentioned above, there are no changes in the solutions for A oo. For A —>■ — oo the 
solution is contained in M: 



sinh((if/i/2T)l/(/i-2)) 



1 2 



sinh((ff/i/2T)(j-l)/(/i-2)) 



j = 3...p-2 



Hp 



For H ^ this implies 



In 



1 



sinh((if/i/2T)(/i - 3)/(/i - 2)) 
sinh(iJp/2T)l/(/i-2) 
sinh((iJ/i/2T)(/i - 3)/{p - 2)) 
sinh(i7/i/2r)l/(/i-2) 



(j - 1) 



2 , j = 3 ... /X - 2, 



1 



1 



/i-2' /i-2 
Now we are ready to find the sum in equation (4.18) 

t(-l)^0)7(^7,^;)=7(l,^)+7(l,^ 



1 



1 1 



7 



7 1,7 +7 1 



1 



jU — 2 yU 

^ r,0 



,(j + iy 



/i-2 



-E7 0, 



(J + 1)^ 



4 



9 



(4.29) 



(4.30) 



- -7 1 



+ 7(0,7) +7(0,-1 +27 



/i-2 /i^ 



1 1 



/i — 2 yU 



14 



TT 



27(l,0) + y-27 

27(1,0) 
27r2 



27 (0, - 



TT 

— + 27 

3 ^ 



Finally, 



Polarization (4.21) is modified 



T 

27rc 



In 



1 + exp{(p 



In 



1 + exp(v9 



1 + exp{-!f: 



l + exp(-(/?+ 

Relations (4.23) are still valid. From equation (4.29) we obtain 

Hfi fi — 3 



Hfi 1 
T ix-2' 



and 



Therefore, 



T 



27rc ^ 
T H 



H 
T 



H - I - fi + 



/i 



27rc T /i - 2 



and finally 



1 



1 



7 



1 1 



— 2' /i 



(4.31) 



(4.32) 



(4.33) 



(4.34) 



(4.35) 



(4.36) 



27rc jj — 2 2ttc tt — 27 

Now we have to compare our results with those of other authors who have presented 
calculations especially for c > 0. To avoid ambiguities we multiply heat capacity and 
susceptibility by the speed of sound Vg, afterwards the result becomes unique, not 
depending on the normalisation of the coupling constant. 
From nil and [0 we can derive 



2c7r 



7 
2\c\tx 



for c > and 



for c < 0. 



TT 



7 



(4.37) 



We have used our method also to obtain the values for the two homogeneous systems 
{s = \ and s = 1). At least, the susceptibility for s = 1 has not been calculated before 
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in the case of negative coupling. (Heat capacity was determined in paper p.) In all 
cases considered we found for the heat capacity the conformal result 

Cv, = ^ (4.38) 

where Cy is the central charge of the Virasoro algebra, which is equal to one for negative 
coupling. 

It is remarkable that formula (4.11) is preserved, because (rewritten for the entropy 
per site ) the factor in front of the sum is always equal 1 / (f svr) (apart from sign) while 
the sum measures the central charge (being equal to 7r^Ct,/3). 

The form of equation (4.21) for the polarization per site can be understood in the 
same way. The factor in front of the logarithms is always T/{2vs'~f), while for different 
signs and spins the logarithms also differ. 

For positive coupling they yield the result 25"/ (tt — 2'~fS')H/T where S" is the larger 
of the two spins (our model has S" = 1), which may be equal. 

Therefore 

^sXic > 0) = ^' (4.39) 
TT — 25'7 

which is consistent with all former results, especially with paper for homogeneous 
chains and with paper |Q for the isotropic limit of alternating chains with S' > S. 
For negative coupling the logarithms always equal H/T leading to 

VsXic <0) = ^ (4.40) 
27 

with no dependence on the spins. This is remarkable, because we remember c„ = 1 in 
the same case. 

5. The isotropic model with alternating spins 

In this section we present some results for the isotropic limit of the model considered 
before, which we will call XXX{^, 1). On one side, there are some peculiarities in the 
limit 7 — > (especially for negative couplings). On the other side, in the sectors with 
different signs of couplings it is possible to obtain several new results yet undiscovered 
for the anisotropic case. 

To begin with we have to define the isotropic limit. The model has been considered 
in paper 0, we have only to specify the normalization of coupling constants to fit with 
our section 2. 

The BAE take the form 

/ Aj + i + A ^ _ il Aj - Afc + 2 



n , j = i...M. (5.1) 
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Instead of (2.3) we define 

CnEn — ( ^ 



E = ciEi + C2E2 -{In-M]H (5.2) 



with 



where 



M 



Ei = -Y.<^j)^ ^ = 1'2 (5.3) 



A2 + ^ 



4 



Taking the hmit of equations (2.4) and (2.5) we see that we have to put 

Ci = hm — and 

7— »o ^ 

Co = hm — . (5.5) 

The TBA (for zero temperature) has been given in 0: 
ei(A) = -27rcip(A) + p*e+(A), 

e2(A) = -27rc2p(A) + p * e+(A) + /i * e+(A) + — . (5.6) 



Here 



p(A) 



2 cosh vrA ' 



= I ^ \^ ^ ^ e'P\ (5.7) 



27r2cosh(p/2) 

As in section 3 we must distinguish between the various regions of signs when the system 
(5.6) is solved. We will follow directly the notation from section 3 above. 

The solutions in the sectors (i) and (iv) are well known. While in (i) the ground 
state is antiferromagnetic (1- and 2-strings) and therefore the limit of the anisotropic 
case, in (iv) the ground state is ferromagnetic and hence different from the anisotropic 
case. This explains why the results of section 3 for c > (with the replacement (5.5)) 
lead to the isotropic values, while they diverge for c < 0. 

Now we wish to investigate in some detail the sector (ii) with ci > 0, C2 < 0. Only 
real roots can be present in the ground state. The TBA (only one equation left) can be 
formulated in two equivalent ways, which we shall need both. 

ei(A) = -27rcip(A) + 27i\c2\h{X) + h * e+(A) + (5.8) 



'^iW = -ciai(A) + |c2|a2(A) - 



Ztt 



H. (5.9) 
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From equation (5.9) one easily determines the region where the solution is ferromagnetic. 
The integral term is always non-negative. Therefore we have ferromagnetic behaviour 
(e(A) > everywhere) if the remaining function of A on the RHS of equation (5.9) is 
strictly positive. That is guaranteed, if it is fulfilled for A = 0. Thus 

- Aci + 2\c2\ + H > (5.10) 

implying 

H^rit = 4:Ci - 2\C2\ (5.11) 

in this region. For vanishing magnetic field ferromagnetism is obtained as long as 

< o < ^- (5.12) 

\C2\ 2 

This is just the isotropic limit of inequality (3.12). For ci = ||c2| there is a phase 
transition to a partially ordered state, the Fermi zone of the 1-strings starts at A = 0. 
The Fermi radius increases and stays finite for i7 7^ moving counter-clockwise towards 
the vertical Ci-axis in the (c2, Ci)-plane. We have strictly proven that there is no point 
where it reaches infinity (for H = 0) unless C2 = 0. One can see that from equation 
(5.9), because for A — >■ 00 h{\) vanishes much slower than p(A). 

Summarizing the facts, (ii) splits into two parts one ferromagnetic (|5.12|) and one 
with partially ordered ground state whose Fermi radius varies from zero to infinity (s. 
figure 2). So one can say that the second coupling C2 works here like an external magnetic 
field rendering the Fermi radius finite as it is for for the homogeneous antiferromagnetic 
models with < H < Hcru- 

Analytical solutions of equations (5.8) or (5.9) can be obtained for large and small 
Fermi radius. We start with the first and consider equation (5.8). It is identical to the 
TBA of XXX(l) model except for the term with C2. We therefore use the technique of 
paper |T^, see also paper |]TT|, recasting that term (after Fourier transformation) as a 
suitable product. We put as usual y{X) = ei(A + b) with ei(6) = y{0) = and use the 
symmetry of ei(A). After Fourier transformation 

/oo 

e^-'f{X)dX (5.13) 
-00 

we write the C2-term on the RHS of equation (5.8) in the form —C{uj)h{uj) with 

C{uj) = -27^|c2|e-^"^ (5.14) 

With the notations of [Tl| the solution in cu-space is given by 

y4uj) = (1 - G4uj))C{uj) + G4uj)Q4uj). (5.15) 

This has to be integrated to obtain the radius b. We have carried out the first integral 
by deforming it along the cut on the negative imaginary axis yielding (for b ^ 1) 

" (1 - G+{uj))C{uj)duj = '^ihjyl. (5.16) 
00 b^ 
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Together with the second part we have the condition 

2vr , |C2|4 




i/-ci2vr^/-e— + ^ = (5.17) 

which determines b = b{H, ci, C2). For its vahdity we have to ensure 

H/ci < 1 and |c2|/ci < 1. (5.18) 
The free energy is calculated via TBA (see e.g. |0): 

— = /o - - J_^p{X) In (1 + e^^W/^) dX-- In (l + e-^/^) rfA. (5.19) 

For vanishing temperature the first integral is proportional to y^{iTi) yielding (in leading 
order) the first term in (5.20). The calculation of the second integral is more involved, 
after one has made use of equation (5.6). We give only the result for the dominant 
terms: 

The term proportional to gives xo = l/(4ci7r^) which in some sense can be 
interpreted as half of the value which we found on the conformal line. This is the 
susceptibility for 

i^<|c2|<ci (5.21) 

because then h does not depend on the magnetic field. The term —H/A describes a 
constant magnetization for \c2\ —>■ 0, which can be found from BAE directly. 

The result for \c2\ -C ^ ci is difficult to interprete, the limit H ^ is not allowed 
in this case. 

Unfortunately we did not succeed in calculating F for T > 0, which meets severe 
difficulties. 

We close the consideration of (ii) calculating the free energy and magnetic 
susceptibility for small Fermi radius, that is close to the line of transition to 
ferromagnetic behaviour. We solve equation (5.9) for H < Hcru with Hcru from above 
(5.11). 

Making an expansion in powers of the Fermi radius b we see that the integral on the 
RHS (except the 5-term) is of power b^. We therefore can easily determine 

ei(A) = er(A) for |A| < 6 (5.22) 

up to terms of power 6^. 
From ei(6) = we have 



Hcrit — H 

16ci + 2c2 



(5.23) 
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and 

ei(A) = (16ci + 2c2)(A2-62). (5.24) 
The free energy per site is given via equation I (3.14) 

and hence 

_ 3 1 1 

This is to be compared with the same value for the usual XXX{^) Heisenberg model 

2 1 1 

Xxxx = ~~/f^~7f T TT (5-27) 

in our normalization of coupling constant. That result is, of course, not the limit C2 — > 
of equation (5.26), because the change in I (3.14) also must be taken into account. 

At the end of this section we shortly comment on the sector (iii). It is treated in 
the same way as above the sector (ii). Equations (5.8) and (5.9) are replaced by 

62(A) = -27rc2p(A) + + { . I + MA)) * 6,- + f , (5.28) 

smh(7rA) \2smh(7rA) J 2 

e+(A) = |ci|a2(A) - C2(ai(A) + a2(A)) - (^5(A) + ?M^L±^j * ^- + 2H. (5.29) 

The critical magnetic field can be read off from equation (5.9): 

8 , , , , 

Hcrit^ ^C2-\C]\. (5.30) 

For vanishing field we have ferromagnetic behaviour as long as 

M > 8 (5.31, 

C2 3 

(compare equation (3.20)). 

The power expansion in h is rather simple while the Wicncr-Hopf calculation for 
large 6 is a little bit more involved. Therefore it is not carried out here. The results will 
qualitatively agree with those from above. 

The phase structure is depicted in figure 2. We have four sectors and two singular 
lines (the positive axes). Three of the sectors show critical behaviour without mass gap; 
one of them is truly antiferromagnetic with infinite Fermi zone. The remaining one is 
ferromagnetic. 
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Figure 2. The phase structure of the XXX{^,1) model. The sector without 
indication of ground state strings is ferromagnetic. Again axes are drawn broken 
except if they coincide with sector borders, where they are drawn sohd. 

6. Conclusions 

Wc have considered the XX Z{^, 1) model with strictly alternating spins in one part of 
its critical region of anisotropy < 7 < 7r/2. We expect a similar but not identical 
behaviour in the other part 7r/2 < 7 < tt because there is no obvious symmetry 
between the two regions. The model contains two parameters, anisotropy and the ratio 
of coupling constants, and shows a rich physical structure. Except for the isotropic 
model (where we have a ferromagnetic region) we found an antiferromagnetic ground 
state and no mass gap. So the model behaves critically, but it is conformally invariant 
only on a line c = c and also in a large sector including this line and having at least one 
negative coupling. Around that line there exist sectors where the ground state does not 
depend on 7, seperated from each other by sectors where it depends crucially on 7. It 
is remarkable that two kinds of sectors are also different with respect to the occurence 
of finite Fermi zones. 

The sectors around the line with equal couplings are well studied now, their ground 
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states and excitations have been established. At the hne c = c we have calculated 
low temperature heat capacity and magnetic susceptibility. Different signs of couplings 
cause very different behaviour, e. g. different central charges (1 or 2) and different 
behaviour of susceptibilities. 

A subsequent paper will deal with finite size corrections in those sectors. We expect 
the standard results in the conformal case while apart from that conformal symmetry 
does not make any prediction. 

The sectors with finite Fermi zones require further treatment, including numerical 
studies. The same applies to heat capacity and susceptibility for different coupling 
constants. 
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Figure captions 

Figure 1. The phase structure of the XXZ{^, 1) model for 7 = n/3. The ground state 
strings are indicated for the four different sectores. The arrows symbohze the decreasing 
Fermi radii of the corresponding strings. Broken lines are coordinate axes. Where axes 
are drawn solid they coincide with sector borders. 



Figure 2. The phase structure of the XXX{^, 1) model. The sector without indication 
of ground state strings is ferromagnetic. Again axes are drawn broken except if they 
coincide with sector borders, where they are drawn solid. 



